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ABSTRACT 



Context. The Crab pulsar emits across a large part of the electromagnetic spectrum. Determining the time delay between the emission at 
different wavelengths will allow to better constrain the site and mechanism of the emission.We have simultaneously observed the Crab Pulsar 
in the optical with S-Cam, an instrument based on Superconducting Tunneling Junctions (STJs) with yus time resolution and at 2 GHz using the 
Nancay radio telescope with an instrument doing coherent dedispersion and able to record giant pulses data. 
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, Aims. We have studied the delay between the radio and optical pulse using simultaneously obtained data therefore reducing possible 
• ■ uncertainties present in previous observations. 

Methods. We determined the arrival times of the (mean) optical and radio pulse and compared them using the tempo2 software package. 

Results. We present the most accurate value for the optical-radio lag of 255±21 ps and suggest the likelihood of a spectral dependence to the 
1 excess optical emission asociated with giant radio pulses. 
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Key words, pulsars: individual (Crab Pulsar, PSR J0534+2200) 

1 . Introduction heights. This should then be reflected in the timing of the pulses 

at different energies. Small differences in pulse alignment allow 



The Crab pulsar was the first pulsar detected as a consequence 
of its occasional bright pulses (now known as Giant Radio 
Pulses) at radio wavelengths (Staelin & Reifenstein 1968; the 

first pulsar ever detected was PSR J1921+2153 Hewish et al. , ,. , , ...... 

i ,„ ,,/T . , . . , spectrum is thus a powerful tool to constrain theories of the 

1968). Later on pulsations have been detected at optical wave- • • • • j i 

lengths (Cocke et al.; |l969l l and throughout the electromagnetic 

spectrum at X-rays and y-rays. One feature which is constant 

over the whole spectrum is the pulsed emission, but the details , , 
,.„. . . .. „ „ . , In recent years, it has become clear that the main and sec- 
diner substantially. E.g. the strength of the interpulse, the pres- , T „„, „„, 

c 7 t , • , , „ , \ ondary pulses of the Crab Pulsar (PSR J0534+2200) are not 

ence of a pre-cursor to the main pulse (at low radio frequen- J r 

. , ... , TT . . . , aligned in time at different wavelengths. X-rays are leading the 

cies) change with wavelength. However, the precise timing of f . J _2 — . 

, a c ac ' radio pulse by reported values of 344±40 us (Rots et al. 2004; 

the pulse over many orders of magnitude of energy imposes se- _„,JT , , ,„ , r . ' 

F . \ . . B . , . i_ \ T7 PvXTE data) and 280±40 us (Kuiper et al. l2003l INTEGRAL 

vere constraints on the emission regions (and mechanism). For , , , , , — T ^ 

in • o ir a- i Jmnci i. . j lU . .u data) and y-rays are leading the radio pulse by 241+29 us 

example Romam&Yadigaroglu 0995) have suggested that the , v J__ t ^ -n.- ™ 

pre-cursor originates at the polar cap, while the pulse and in- 
terpulse originate in the outer gap in the magnetosphere, with 
higher energy pulses being generated at significantly greater 



to study the emission regions at relatively small scales (about 
10 km per 0.001 phase in case of the Crab pulsar), and precise 
timing of pulsar light curves throughout the electromagnetic 



spatial distribution of various emission regions: a time delay 
most naturally implies that the emission regions differ in posi- 
tion. 



(Kuiper et al. 2003 1 EGRET data. The uncertainty in this value 
does not include the EGRET absolute timing uncertainty of 
better than 100 /-is). At optical wavelengths, the observations 
presented a less coherent picture. Sanwal (1999) has reported 
Send offprint requests to: T. Oosterbroek a time shift of 140 fis (optical leading the radio). The uncer- 



2 



T. Oosterbroek et al.: Simultaneous Absolute Timing of the Crab Pulsar at Radio and Optical Wavelengths 



Table 1. Log of optical and radio observations. All times are in 
UT. The first two optical observations were performed with a 
different instrumental setting resulting in different instrumental 
delays (see text). 



Observation 


exposure (s) 


Start time (yyyy-mm-dd hh:mm:ss) 


0001 


179 


2007-09-16 04:05:49 


0002 


900 


2007-09-16 04:09:05 


0003 


900 


2007-09-16 04:25:58 


0004 


900 


2007-09-16 04:41:03 


0005 


900 


2007-09-16 04:56:09 


0006 


900 


2007-09-16 05:11:14 


0007 


900 


2007-09-16 05:26:19 


0008 


900 


2007-09-16 05:41:25 


0009 


373 


2007-09-16 05:56:30 


0013 


1800 


2007-09-17 04:10:47 


0014 


1800 


2007-09-17 04:41:42 


0015 


1800 


2007-09-17 05:11:48 


0016 


1100 


2007-09-17 05:41:53 


Radio 1 


3420 


2007-09-15 05:17:10 


Radio 2 


3930 


2007-09-16 05:12:50 


Radio 3 


3960 


2007-09-17 05:08:40 


Radio 4 


3890 


2007-09-18 05:04:50 



tainty in this value is 20 [is in the determination of the optical 
peak and 75 jus in the radio ephemeris. Shearer et al. ([2003 ) 
have reported a lead of 100+20 //s for simultaneous optical and 
radio observations of giant radio pulses. Golden et al. ([2000) 
have reported that the optical pulse trails the radio pulse by 
~80±60 /is. Romani et al. ([2001 ) conclude that the radio and 
optical peaks are coincident to better than 30 /is, but their er- 
ror excludes the uncertainty of the radio ephemeris (150 yus). 
Oosterbroek et al. (2006, hereafter O06) have studied this issue 
in detail and found an optical lead of 273±100 //s. In this paper 
we will use simultaneous optical and radio observations at 2 
GHz to further improve the accuracy of this value. 

We also study the possible changes in the peak emission co- 
incident with Giant Radio Pulses (hereafter GP) as previously 
reported by Shearer et al. (120031) . 

2. Observations 

Optical observations were obtained with S-Cam3 (Martin et al. 
2006 ) mounted on the ES A Optical Ground Station (OGS) tele- 
scope on Tenerife on September 16 and 17 2007. For test pur- 
poses the first two observations were obtained in the so-called 
FAST mode of the instrument. This mode only differs in the 
digital filter used to obtain the pulse shape of the photons. In S- 
Cam3 each photon gives rise to a bi-polar signal in the detector 
electronics chain, the zero-crossing of this signal is used to time 
tag the photon. The instrumental delays in both modes were 
determined by triggering an LED using a GPS and registering 
the assigned time of the observed photons (see also O06). The 
instrumental delay of S-Cam3 amounts to 39.7±0.5 fj.s in the 
instrumental mode which was used for all but the first two ob- 
servations, and 14.0±0.5 yus for the first two (i.e. "FAST mode") 
observations. Note that these values are different from those re- 
ported in O06, since a physically different detector array and 
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Fig. 1. The pulse profile template used for obtaining the TOAs. 
The profile shown is the modestly smoothed profile (see text) 
which at this scale is indistinguishable from the unprocessed 
data. In the zoom-in the difference between the data points (cir- 
cles with error bars) and the smoothed data (line) is shown for 
the peak where the differences are largest. At the bottom of 
the main plot the difference between the smoothed and non- 
smoothed data is shown. 

different instrument modes were used. Observations 0009 and 
0016 have a shorter duration than the typical exposure time 
for the night since they were interrupted because of rising sky 
background (twilight). Some thin cirrus was present at the end 
of the first night (starting at around 5:30 UT), while the see- 
ing was somewhat better during the first night compared to the 
second night. 

Radio observations at 2 GHz were performed using the 
Nangay radio telescope on September 15-18 2007. The obser- 
vations were done over a bandwidth of 64 MHz (from 2016- 
2080 MHz) in 16 channels. Data were coherently dedispersed 
over 4 MHz bandwidth channels (providing 250 ns resolution) 
using a transfer function in the complex Fourier domain of 
the recorded voltages (Cognard & Theureau 120061 Demorest 
2007). Using a polynomial representation for the expected evo- 
lution of the apparent rotational phase of the pulsar, data were 
folded to produce a profile every minute. A daily profile is built 
by integration of the 1 -minute profiles obtained over the typi- 
cal 1-hour duration of an observation. During data processing, 
the instrumentation is able to detect and store shortlived events 
such as giant pulses known to occur on this pulsar. Radio ob- 
servations at this high frequency allow to minimize the effect 
of uncertainties in the Dispersion Measure (DM). For a log of 
the optical and radio observations we refer to table Q] 

3. Analysis and Results 

3.1. Optical TOAs 

The analysis of the optical data consisted of searching for 
the best local period by using a profile-folding technique for 
each S-Cam observation. Since these observations are rela- 
tively short (up to 1800s) we did not take into account the pe- 
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riod derivative nor did we perform any barycentering on the 
time stamps. The data was then folded on the best local pe- 
riod, using the midpoint of the observation as the epoch. This 
resulted in 13 profiles, one for each observation. A template 
profile was constructed. This template was created by apply- 
ing a barycentric correction to the time stamps of the pho- 
tons and folding on phase (i.e. calculating the phase of every 
photon using an epoch and period and summing to obtain a 
phase histogram) using the best period determined by folding 
the data on a range of trial periods and maximizing the \ 2 with 
respect to the mean. All data were used, except the two first 
observation, since they are obtained in a different instrumen- 
tal mode. A cross-correlation between every observation and 
the template profile was performed and the delay/lead was de- 
termined by fitting a Gaussian around the peak of the cross- 
correlation curve. The time of arrival of each average pulse 
profile was then assigned to the midpoint of each observation. 
These TOA's were corrected for the instrumental delay of S- 
Cam (14.0+0.5 [is for the first two observations, 39.7+0.5 [is 
for the remainder, see Sect. [2]). Using the cross-correlation tech- 
nique the most significant contribution to the uncertainty in the 
TOA is the uncertainty in the determination of the peak of the 
template profile, which amounts to 18 [is. This uncertainty was 
determined by fitting the data points with a Gaussian in a 0.01 
wide phase range around the peak. This phase range was deter- 
mined iteratively and choosen symmetrically around the peak 
(see also O06). The uncertainty represents the 1 cr uncertainty 
on the peak position (Ax 2 = 1). This uncertainty is larger than 
the uncertainty introduced by the cross-correlation since for 
the cross-correlation the whole profile is used. The individual 
uncertainties from the cross-correlation are an estimated 1 .4— 
3.3 [is. The value of 18 [is can be considered as a systematical 
error on all the indivual optical points, since they are all deter- 
mined using the same template profile. 

In order to test whether some of the scatter in arrival times 
was caused by the presence of noise in the template profile we 
constructed a template profile using all data (except the first two 
observations) which was modestly smoothed by a Savitzky- 
Golay filter (see Fig.[T|). The filter was constructed using 4 for- 
ward and 4 backward data points and degree 4 (see Press et al. 
1992). This almost did not affect the results: arrival times did 
not change by more than 4.3 [is and the change was typically 
less than 2 jus in a seemingly random way. This implies that the 
obtained time shift does not depend on some modest smoothing 
of the template pulse profile. 

3.2. Radio TOAs 

The radio time of arrivals (TOAs) were derived for each day 
of observations yielding 4 TOAs. Each TOA was obtained by 
doing a cross-correlation of each daily profile with a mean pro- 
file obtained by integration over many months of observation. 
This cross-correlation gives a time and a uncertainty ranging 
from 3 to 8 [is. The relatively low signal-to-noise ratio (around 
20) for the daily profiles explains the limited accuracy of the 
TOA determinations. The reference point in the mean profile 
was chosen to be the phase bin with highest radio flux. Given 
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Fig. 2. Example of a giant pulse detected on Sept 17th with 
the Nancay radiotelescope. Total intensity is plotted for every 
4 MHz channels with a resolution of 250 ns. Duration of this 
giant pulse is of the order of 2 [is. 

the signal-to-noise ratio and the width of the phase bins, we es- 
timate that this introduces an uncertainty of 8 [is which affects 
the ensemble of radio TOA in an identical way (and does not 
disappear when averaged over multiple TOAs). 

The determination of the radio TOAs included the Nancay 
instrumental delay. This delay was estimated from several ob- 
servations of the millisecond pulsar J 1824-2452 done with the 
Green Bank and the Nangay radiotelescopes at the same epoch. 
A careful analysis of the mean profiles obtained with the two 
similar instrumentations installed at Green Bank and Nancay 
(both Astronomy Signal Processor -ASP-, Demorest, 2007) 
give a relative offset of 7.0+0.4 [is. Knowing accurately the 
Green Bank instrumental delay (which amounts to 15.74 [is; 
Demorest, 2008), we can derive the Nangay instrumental delay 
as 8.7+0.4 [is. 

3.3. Giant radio pulses and DM determination 

Accurate determination of DM variations is usually done from 
TOAs obtained at two different frequencies, as much sepa- 
rated as possible in order to maximize the effect of the dif- 
ferential dispersion delay to measure. However, as shown by 
Ramachandran et al. (2006), highly frequency-dependent mul- 
tipath propagation in the turbulent interstellar medium pro- 



T. Oosterbroek et al.: Simultaneous Absolute Timing of the Crab Pulsar at Radio and Optical Wavelengths 



4 



CO 



E 



54357 54358 54359 54360 54361 54362 

Date (MJD) 

Fig. 3. The DM value determined from the Nangay GP data. 
Displayed are the nightly mean values. 

duces huge variations of the space (cigar shape like) within 
which the radio waves senses the electron density. TOAs ob- 
tained at very different frequencies did not see exactly the same 
electron density and should not be used to derive a highly re- 
liable Dispersion Measure value. Using frequencies as close as 
possible tends to remove this effect but then TOAs determina- 
tion must be much better to compensate. As Giant Pulses TOAs 
are characterized by very short bursts of radio emissions suit- 
able for a very good TOA determination, we should be able to 
use them even on limited frequency span to determine unbiased 
DM variations. For example, a differential delay between 2016 
and 2080 MHz known at 0.5 /is gives roughly the same DM 
uncertainty as a 10 //s accuracy between 1400 and 2050 MHz. 

We have determined the DM value from the Nangay data 
alone in the following way: each GP (an example is shown in 
Fig. [2]) which was detected in more than 14 channels (among 
the 16 available, using a 15 sigma threshold for each channel) 
was dedispersed with different DM values and integrated in 
frequency to produce a 56 to 64 MHz total bandwidth profile. 
Based on the fact that a good DM value will tend to properly 
align the individual GPs and thus maximize the amplitude of 
the integrated profile, this maximum amplitude was kept for 
each trial DM and a gaussian fit was done on those maximum 
amplitudes to find the best DM. An average of the numerous 
DM obtained for a given observation (one per single GP event) 
provides us with a nightly DM value (see Fig.0. Uncertainties 
are estimated from the rms of the different DM values obtained 
each night. The overall average value over the 69 DM values 
obtained over 4 nights is 56.7740 pc.cm -3 and again an uncer- 
tainty of 0.005 is estimated (from the rms over the 69 individual 
values). 

3.4. Optical-radio lag 

We determined an improved solution to the v and v starting 
from the Jodrell Bank ephemeris to the TOAs of the mean pulse 
using the Nangay radio data. While the fit to the 4 data points is 
not very constraining, it establishes a reliable baseline against 



which the optical TOAs can be compared. The obtained values 
are v = 29.755351174(9)s" 1 , andv = -3.727(6) x lO^V 2 (for 
a fixed value of v of 3.72 x 10~ 22 s~ 3 ). The dispersion measure 
was fixed at the value reported above of 56.7740 pc cm -3 . 

The OGS geocentric coordinates used in the reduction 
(within tempo2) are (X,Y,Z) = (5390280.91, -1597890.59, 
3007083.48). These coordinates were obtained from the read- 
ing of two different GPS units (see O06). The obtained optical 
TOAs were then used in the tempo2 analysis package (Hobbs et 
al.|2006), together with the radio TOAs. The radio and optical 
data points were fit together allowing for an offset between the 
radio and optical. Since a comparison of arrival times is made 
at (almost) the same time, this procedure is not substantially 
affected by e.g. uncertainties in the pulsar position (proper mo- 
tion) or uncertainties in e.g. frequency and derivatives. 

In Fig.@]the results of this analysis are displayed: an offset 
of 255+4 yizs between the radio peak and the optical peak. The 
uncertainty in this value is purely the statistical uncertainty in 
the best-fit value for an offset between the radio and the optical 
measurements. This value does not take into account uncer- 
tainties in the DM values, which could be of the same order. A 
difference in the DM value of 0.005 pc crrT 3 (more than a typ- 
ical month-to-month variation, see Fig. [5j will result in a time 
difference at 2 GHz of ~6 [is. We therefore assume an uncer- 
tainty of 6 [is as a result of uncertainties in the DM value. This 
uncertainty in time is smaller than quoted by other authors (e.g. 
Lyne et al. 119931 1. since our radio observations have been ob- 
tained at a higher frequency which is substantially less affected 
by timing delays due to scattering. 

Adding all of the above uncertainties (18 yt/s from the op- 
tical template, 8 [is from the radio template, 6 [is from DM 
uncertainty, 4 [is from the fit to the offset between the radio and 
optical TOAs) linearly we obtain a total uncertainty of 36 [is. 
This is likely to be an overestimate of the total uncertainty since 
all individual deviations have to be in the same direction to end 
up at this value. Adding all uncertainties in quadrature yields a 
final uncertainty of 21 /is. 

From Fig. [4] we conclude that the variation in the optical- 
radio delay is substantially larger than the individual uncertain- 
ties on the TOAs (which are between 1.4 and 2.1 [is, \<r). The 
average value of the delays obtained during the first and second 
night differ by 6 [is (only 1 .4 [is when one excludes the deviant 
point in the second night). The observed scatter during the first 
night and the possible difference between the two nights, sug- 
gests that small variations (less than about 10 [is) exist. For this 
relative comparison the uncertainties introduced by determin- 
ing phase zero on both the radio and optical templates (which 
each amount to 8 [is) should not be considered, since they af- 
fect all data points in the same way. We consider it unlikely that 
this night-to-night variation is solely caused by variations in the 
DM, since a variation of ~1 times that of typical monthly varia- 
tions is needed in 1 day. Variations in the peak position on time 
scales of days and shorter are studied in more detail by Karpov 
et al. (2007), who find, during one observation, variations at the 
~150 jt/s-level. 
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Fig. 4. Comparison of radio and optical results. At the bottom 
a zoom-in on the optical data points only is presented. Plotted 
error bars represent 1 cr uncertainties. 
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Fig. 5. The recent (past 1000 days) history of the dispersion 
measure towards the Crab pulsar (data taken from Lyne et al. 
2007 ). The data points are monthly values. The time of our ob- 
servations is indicated by an arrow. The difference between the 
Nangay DM value and the value plotted here is less than 2 <x 
(offset is 0.008 while uncertainty is 0.005 for both determina- 
tions). As the frequency spans used for these two derivations 
are not the same (610-1400 MHz for Jodrell and 2018-2080 
MHz for Nan§ay), the difference is acceptable (see text). 




3.5. Giant radio pulses and optical 

In total 927 Giant radio pulses were detected during the radio 
observations. The rate of GPs was much higher during the sec- 
ond night where we have optical data (on average 1 per 6.6 s) 
than the first night (1 per 29.5 s). 

Following the method by Shearer et al. (2003 ) we have se- 
lected the optical data coincident with Giant Pulses (GPs) de- 
tected in the radio observations. The barycentric arrival times 
of the GPs were computed using tempo2, while the event time 
tags were corrected using our barycentering code. We checked 
that these corrections are consistent to better than 15 yus and 
that therefore this procedure provides the required accuracy. 
We find, in total, 603 GPs coincident with optical data. We 
made sure that these GPs occurred near the main peak of the 
profile and selected data half a phase range before and half a 
phase range after the main peak (in the radio). The data was 
then folded to obtain a pulse profile coincident with radio GPs. 
A reference profile was also constructed from data spanning 
80 cycles around (but not including) this cycle. This was done 



Fig. 6. A comparison of the optical pulse profile obtained co- 
incident with a Giant radio Pulse in red and a reference profile 
in green (see text). The difference between the two is also dis- 
played. 



since changes in sky background and airmass force us to select 
data close to the cycle coincident with the GP. Several different 
ranges were tried and all gave similar results, with obviously 
better statistics as the number was increased. We settled for 80 
cycles, since this gives a very good reference profile and yet it 
still only spans 2.7 s of data. We ensured that no optical data 
was missing during these periods (e.g. at the end of an observa- 
tion). We present the obtained results in Fig. [6] As can be seen 
from this figure no clear excess is present. Depending on the 
choice of phase range this excess amounts to about ~ 2.1cr. We 
view this a confirmation of the result of Shearer et al. (|2003 ) 
and shows that their dataset has a much better statistical quality 
given the longer exposure times and larger telescope aperture. 
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Fig. 7. The delay between photons in different energy bands. 
Wavelengths indicate the center of the three bands (not the av- 
erage energy of the photons). 
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Fig. 8. The radio delay with respect to other frequencies. The 
S-Cam data point is presented in this paper, the RXTE value 
is from Rots et al. d2(J04l the INTEGRAL and EGRET points 
are taken from Kuiper et al. (2003). We have added a 100 //s 
error bar to the EGRET value, because of uncertainties in the 
EGRET absolute time calibration (see Kuiper et al. 2003 ). 



3.6. Wavelength dependence of the delays 

In order to investigate whether we are able to detect a delay be- 
tween the pulse profiles obtained at different wavelengths we 
divided the photons obtained during the observations described 
above in three adjacent bands, each containing approximately 
the same number of photons, using the intrinsic energy resolu- 
tion of S-Cam. The wavelength boundaries of the three bands 
were: 3920-5490 A, 5490-6330 A, 6330-8230 A. We then cre- 
ated the folded profiles per observation and did the crosscor- 
relation with the smoothed profile (which was obtained using 
all photons) as described above. The results are displayed in 
Fig. [7] While the quality of the data does not seem to justify 
drawing strong conclusions, a linear fit to the datapoints re- 
veals that the slope in the delay is -5.9+1.9 /is/1000 A, which 
seems to indicate that the blue photons arrive, averaged over 
phase, slightly earlier than the red photons. A similar result is 
obtained when only the main peak part of the pulse profile is 
used for this analysis. It is important to realize that this is only 
a statement about the phase-averaged arrival times of photons: 
the full width at half-maximum (FWHM) of the pulse peak is 
known to decrease with decreasing wavelength (Percival et al. 
[19931 Eikenberry et al. [1997l Sollerman et al. [20D0b . This trend 
is also evident in our data. Together with the assymmetry in the 
profile the decrease in FWHM with decreasing wavelength will 
lead to an average time shift. Close inspection of our profiles in 
colour suggests that the profile rises slightly earlier in the blue 
band, which could yield the result present above. Furthermore 
the separation between the main and secondary peak is smaller 
in the UV than in the optical (Percival et al. 119931 ). In essence 
this shows the limitation of using the cross-correlation method: 
it is only possible to obtain a proper time shift if the profile and 
template are identical in shape. Any other method (e.g. fitting 
around the peak) will also be affected by the assymetry of the 
profile and will give similar results. 



4. Discussion and Conclusions 

We have performed simultaneous observations of the Crab pul- 
sar in the optical (3500-8200 A) and radio (2 GHz). The re- 
sults provide the current best estimate of the optical lead (with 
respect to the radio peak) of 255 ±21 fi&. The two major contri- 
butions to the uncertainty are the determination (or definition) 
of the peak in the optical, since the profile is asymmetric, and 
uncertainties in the DM value which translate into significant 
uncertainty in the correction of the radio data. Additionally a 
true jitter in the optical peak on different time scales might 
take place with an amplitude of 5-10 //s. However, our value 
is an order of magnitude more precise than determined by O06 
and fully consistent with the value reported there. The close 
match with their value (273+100 //s) must be fortuitious given 
the large uncertainties on the values, and can be considered a 
statistical coincidence. A comparison with previous results can 
be be found in O06, however due to the increased precision 
we can now conclude that the result of Shearer et al. (|2003), 
an optical lead of 100+20 fi&, is not consistent with our result. 
Together with the result of Golden et al. (2000), an optical de- 
lay of ~80±60 jus, which was already inconsistent with the re- 
sults presented in O06, this raises the possibility of a long-term 
change in the radio-optical difference. 

From a plot of the radio delay throughout the electromag- 
netic spectrum (see Fig. [8]) it appears that this delay is ap- 
proximately constant over more than 7 orders of magnitude. 
We note that the non-optical results have been derived from a 
comparison with, non-simultaneous Jodrell Bank observations, 
which results in larger error bars. This near-constancy of the 
delay with respect to the radio profile implies a similar emis- 
sion mechanism and emission location from the optical to y- 
rays. This is indeed the case in modern models by e.g. Takata 
& Chang (120071 

Our results for the stability of the pulse profile (differ- 
ences of up to ~10 jjs on a time scale of a day and even 
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smaller variations on a time scale of less than a day; see Fig. [4]> 
are not consistent with the results from Karpov et al. (2007) 
who present variations of up to ~150 jus in a fraction of a 
day. However, this variation is only present in a small part 
of their data set (December 2005 and January 2006), while 
the remainder of their data (obtained in December 1994 and 
1999, and November 2003) do not show substantial variations. 
This suggest that the behaviour observed in December-January 
2005/2006 is an anomaly and that night-to-night variations in 
the Crab optical peak positions amounts to at most ~10 yus. 

Our value for the time delay between red and blue pho- 
tons of -5.9+1.9 jus/1000 A can be compared to the results of 
Sanwal ( 1999) who found a delay between the R and U band 
of -10+4 fis and Golden et al. (120001) who found that the peaks 
in B and V are coincident within 10 fis. Taken together these 
three measurements indicate a small delay between blue and 
red photons. However changes in the FWHM of the pulse pro- 
file and the assymmetry of the profile and possible differences 
in the separation between the main and secondary peak make 
an interpretation complex. 

Current models for high-energy pulsar emission (e.g. Dyks 
& Rudak 2003), describe the peaks as caustics: special relativ- 
ity effects (aberration of emission directions and time of flight 
delays due to the finite speed of light) cause photons emitted 
at different altitudes to pile up at the same pulse phase. It is 
difficult if not contrived to imagine how caustics originating 
from the trailing edge extending along the full length of the 
open field lines within the magnetosphere could produce phase 
bunching yielding a common main peak cusp spanning 7 orders 
of magnitude in photon energy from incoherent synchrotron 
emission, and less than 1 % in phase difference from coherent 
radio emitting sources, spanning in this case some 12 orders in 
photon energy. As such arguments are usually presented in the 
context of outer gap or slot gap model formalisms, these new 
results argue strongly for a localised source for the observed 
multiwavelength emission from the Crab pulsar. 

That the radio and higher energy emission process are in 
someway linked is not in doubt, following on from the optical- 
GP association originally shown by Shearer et al. (2003). 
Our confirmation of an excess, albeit at low significance fur- 
ther strengthens the connection between these energy bands. 
Variations in the optical emission correlated with GPs have 
been argued to be a consequence of local density enhancements 
to the plasma stream (Shearer et al. 2003). Recent studies of 
GPs from the Crab have indicated that these can and do orig- 
inate as short-lived (~ ns), narrowband nanoshots, which have 
been argued to have their origin in strong, highly localised 
plasma turbulence (Hankins & Eilek, 2007). That there is no 
change in phase or shape of the Crab's optical light curve mor- 
phology during a GP associated event, just flux, requires an ex- 
tended region of plasma turbulence to produce the same phase- 
bunching effect following the caustic model ansatz, which is 
clearly implausible. The measured time delay across the opti- 
cal bands represents < 3.0 X 10" 4 of phase, with some evidence 
of a jitter in the optical peak across several time scales occur- 
ing with a similar variation in terms of phase. Again, these ef- 
fects are not associated with any measurable change in incident 
flux. Assuming a constant emission source, localised geomet- 



ric constraints represent the most likely cause of such observa- 
tional effects. Together these interpretations merit some revi- 
sion of the fundamental approaches to modeling the Crab pul- 
sar's emission processes, and there may be some use in trying 
to geometrically constrain emission regions within the magne- 
tosphere first, and then address the local physical conditions 
that are consistent with the data we have available to us - rather 
than the other way round. 

Finally the anticipated detection of excess optical emission 
associated with GPs extrapolating from the work of Shearer et 
al. (2003) was expected to be < lcr prior to our observations, so 
its marginal detection is somewhat surprising. We note that our 
passband (~3500 — 8000 A) is quite different from the pass- 
band (6000-7500 A) of Shearer et al. d2003l and might result 
in different results depending on the spectrum of the "excess" 
optical emission. If our observed excess had been significant 
we could have constructed the optical spectrum obtained when 
a GP occurs and compare it with the reference spectrum when 
no GP is present. This is possible by using the intrinsic en- 
ergy resolution of S-Cam (each photon is tagged with a time 
stamp, but also its energy, Martin et al. 120061 ). Such an ap- 
proach could reveal the spectrum of the excess emission and 
consequently would impose important further constraints on 
the emission mechanism. We regard the acquisition of deeper 
coincident radio-optical data as critically important to confirm 
the spectral dependence of such an excess. 
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